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Abstract 



We describe a first attempt to calculate scalar 2-loop box-functions with arbitrary internal 
masses, applying a novel method proposed in Q. Four of the eight integrals are accessible 
to integration by means of the residue theorem, leaving a rational function in the remaining 
variables. The result of the procedure is a three- or sometimes two-dimensional integral 
representation over a finite volume that can be further evaluated using numerical methods. 

^ , 

The study of higher-loop quantum corrections in QFTs is accompanied by rapidly increasing 
computational challenges. While in pure QED e.g. the presence of only one mass scale allows 
relative high loop orders to be evaluated , the feasibility of loop calculations in the general case 
with arbitrary internal masses has up to now been cut off at the two-loop level. 

In recent years much work has been done on two-loop integrals with arbitrary internal masses 
of the self-energy- and vertex-type. These functions could be reduced to double integrals with 
manageable numerical behaviour J|, [| ||, U ( c -f- @> Hi ^ or an °ther approach at these diagrams). 
Not much seems to be known, however, about the different two-loop corrections to scattering- 
amplitudes (Fig. [l]) involving arbitrary-mass propagators which will become of some importance for 
probing new physics in W + W~ -processes or in a Muon-collider, for example. Some of these topolo- 
gies factorize into products of one-loop functions, i.e. no propagator shares both loop-momenta. 
Others are one-loop functions containing other one-loop functions as subtopologies (Fig. |l|, third 
column). Rather than analyzing all of these functions, we focus our attention on five four-point 
topologies. They are all obtained by shrinking propagators in the two basic four-point topologies 
I X and II. Two of them are of crossed type where two propagators share both loop-momenta 
(first column) and three of planar type (second column). There is no striking difference to the 
graphs in the third column. For example the first graph in the first row, third column 5 relates 
to the graph in the second row, second column H if we shrink the appropriate propagator. 

The importance of the four point function derives not only from the fact that it describes 
the first scattering amplitude (2 — > 2), but also from the fact that the multiparticle scattering 
amplitudes, involving n-point functions, n > 5, are closely related to the four point functions ||. 
We will see below how this relation emerges in our method at the two-loop level. 
The idea behind our method is to perform four of the integrals with Cauchy's residue theorem 
and evaluate at least one of the remaining ones analytically, leaving a two- or three-dimensional 
representation for numerical evaluation. 

This article is organized as follows: In the first two sections we describe how one can generally 
perform four integrations with the residue theorem alone and analyze the structure of the result. 
In the remaining sections we will complete the calculation for a graph of type /V in a special 
kinematic case, always keeping in mind the more general graphs, notably I and IX as we go 
along. In section || we look at some of the results and discuss what possible problems we need to 
be concerned with in some limiting kinematical regimes. 
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Figure 1: Possible 2-loop box-topologies. Starting from class I, where all topologies are constructed 
by appending 4 external legs to the 2-loop master-topology, one can construct all other box- 
topologies by successive cancellation of internal propagators. 

1 Integrating the middle variables 

The functions under consideration can generally be written as 

v = I d*k [ dH - — 1 , (l) 

with five to seven inverse scalar propagators of the form Pz, mi = (l+p) 2 — mf +irj, I denoting a loop- 
momentum and p some combination of external 4-momenta. Note that we do not attach an index 
% to the r\ from the Feynman prescription since we will choose them all to be equal. V depends on 
six independent kinematic variables, a possible choice being the Mandelstam-variables s, t and u 
together with the masses rrii of the four external particles and the condition s + t + u = J^i m i ■ 
We will, however, need to choose an explicit Lorentz frame for our purposes. 

When trying to perform one of the integrals using the residue theorem, it is very attractive to 
first linearize all the propagators Pi with respect to the corresponding variable since this makes the 
detection of poles particularly simple and does not introduce unnecessary square roots inhibiting 
further integrations. Due to the signature of the Minkowski metric this linearization can easily be 
done in one pair of variables. Choosing l\ and k\ for linearization and applying the shift 

h — >lo + h, k a — >fco + fci, (2) 
a propagator Pi undergoes the transformation 

{I + p) 2 -m 2 + iT]= (l +Paf - (h + Pi) 2 - {h +P2? - ih + Psf - m 2 + if] 

— ► (lo + Pof + 2h{h+Pn-pi)-p\- (h + P2) 2 - (h + P3) 2 - m 2 + irj. (3) 

This shift can safely be done since the functions defined by topologies la), II a) , and lb) - III b) 
converge absolutely as can be seen by counting powers. 

Closing the contour in the upper half-plane and using the residue theorem to carry out the l\- 
and fci-integrations, one obtains constraints for the lo- and fco-integrations. This appearance of 
constraints is due to the position of poles in the complex l\- and fci-planes either in the upper or 
the lower half and hence either contributing as a residue or not. These constraints affect the Iq- 
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and /co-integrations only because no other loop-variables appear in the Zi-linear term in (||). We 
will return to these constraints in the next section. 

Noticing that the integrations over the variables h, h, and are still unbounded, suggests 
solving two of them with the residue theorem again. In Q] it was shown how the linearization 
necessary for simplifying this task can be carried out in what we may call the middle variables: 
ki, l\, ki and h- We will briefly repeat the argument. 

As a consequence of conservation of four-momentum, the external legs of any four-point func- 
tion span a 3-dimensional subspace of momentum-space — the so-called parallel space. Using 
Lorentz-invariance, its complement — the orthogonal space — can always be chosen to be paral- 
lel to the 3-axis. The 3-components of the loop- momenta do not mix with any external momenta 
then: 



Pi = l 3 + (something real) + it], 
Pi+k = (h + kz) 2 + (something real) + ir], 
Pk = k 2 + (something real) + ir]. 

Hence the poles of the integrand in the complex I3, k% and (£3 + fc3)-planes respectively are all 
located in the first and third quadrant. 

At this stage one can make contact to the case of a general n-point function. For n > 5 
the 3-components do mix with external momenta. But the very fact that n > 5 ensures that 
after undertaking appropriate partial fraction in the propagators the convergence of all integrals is 
sufficiently strong so that termwise linear shifts are allowed which erase the appearance of external 
momenta in the 3-components of all quadratic propagators. 

The observation that the poles of the integrand are in the first and third quadrant (together 
with the fact that the integrand falls off sufficiently rapidly at large I3, k^ and h + k^) suggests 
that we should try to rotate clockwise by 7r/2 and effectively change the metrics from the usual 
Minkowski metric (+, — , — , — ) to (+, — , — , +) in order to be able to linearize the integrand in I2 
and fc 2 . The mixed propagator Pi+ k , however, seems to spoil this project since its roots in the 
complex £3- or ^-plane alone are not bound to the first and third quadrant. In order to do the 
rotation, we have to treat I3, fc 3 and I3 + fc 3 on the same footing. Due to the integrand's symmetry 
we can restrict our attention to the first quadrant in the fc3-Z3-plane: 
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Jo ■ i \Pl(l 2 3)Pl +k ((l3 + k3) 2 )P k {k 2 3 ) 
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P l (ll)P l+k {(l3-k3) 2 )P k {kl)--- ) 

Now we reparametrize this quadrant by substituting I3 — > uv 2 and k 2 — > (1 — u) v' 2 



1 f 1 du f + °° , ( 1 
v dv 



2 Jo y/u(l-u) Jo V p ^( ^t ' ^;2 ) p ^+fc((^/" + ^/ ^31I ) v 2 ) Pk({l ~ u)v 2 ) 

1 



P t (UV 2 ) P l+k ((y^ - VT^H) 2 V 2 ) P k ((l - U) V 2 ) 



Here, u, (1 — u) and (\/u — y/1 — u) are all positive so the propagators have their poles in the first 
and third quadrant of the complex w-plane. This allows us to close the contour of u-integration 
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around the fourth quadrant, with the Jacobian picking up a sign: 

1 f 1 du f + °° , ( 1 
— v dv 



2 J y/u{l-u) Jo \ Pi(-uv 2 )P l+k (-{^I+VT^ll) 2 v 2 )P k (~(l~u)v 2 ) 



P^-uv 2 ) P l+k (-(^-VT^) 2 v 2 ) P k (-(1 - u) v 2 ) 



Inverting the above transformations we are left with the identity 



+oo r+oc -| 

dl-i / dk 



Pi(l 2 3 )Pi + k{(h + h) 2 )Pk(k 2 )--- 

^ dh J ^ dk 3 p ; (^2)p ;+fc (_ (/3 + fc3) 2 ) p fc (_ fc 2 ) ...- ( 4 ) 

Note that this flip in the sign of metric is not related to the standard Wick rotation where an 
appropriate analytic continuation has to be applied at the end of calculation in order to obtain 
the Greens function for arbitrary exterior momenta. The flip Q is a simple analytical formula 
not for our case only but for any loop-component belonging to orthogonal space. In particular, 
we are not allowed to set the imaginary part in the propagators to zero. 

Now we are in a position to complete the linearization of our propagators in the middle variables 
ki, l\, k 2 and l 2 by applying the shifts 

^3 — >^3 + ^2, k 3 — >k 3 + k 2 (5) 

in addition to (Q). 

Next we interchange the order of integration in order to apply the residue theorem to the four 
middle variables first. Recall that we are allowed to interchange the order of integration for our 
graphs because the integral over the modulus of their integrand exist. 

As mentioned above, the sign of the linear coefficient of the variable being integrated (ii in (|^)) 
determines whether the pole is inside or outside the contour. We therefore obtain a sum of residues 
each having a Heaviside function constraining the domain of integration of the edge variables ko , 
lo, k 3 and I3 arising at each integration of the middle variables. 

2 Constraints 

At each integration, we expect to obtain one residue for each propagator containing the integration- 
variable. The potential proliferation of terms is, however, drastically reduced by two different kinds 
of relations holding among them as we show next. 

The first relation is a well-known corollary to the residue theorem (see e.g. [p^| ): 

Lemma 1 The sum of the residues of a rational function (including a possible residue at infinity) 
is zero. 

Therefore we can express one of the remaining residues after each integration by the sum of all the 
other ones. In our case we apply this to a product of inverse propagators linear in an integration- 
variable I: Pi t i = a>i + Pil + irj, with 77 > 0, n > 1 and containing the remaining mass- and 
momentum-terms. We are allowed to drop one term and only keep its 6*-function: 

J_ f +OC d l 1 - fl(-li) - g(-ft) 
2mJ_ 00 Pi A Pi, 2 PiMmPi 



2ttiJ_ 00 Pi A Pi, 2 Pi,3 PiMm p tMm Pi PiAm Pi,z\m 



-A) - fl(-AQ 



(6) 
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where iW denotes the zero of Pi t i. Since this collecting can be done in four consecutive integrations, 
one can reduce the number of terms for the planar box-function I from 108 to 36. 

The second relation holding among the terms is a consequence of the consecutive integration 
in l\ and 1% and can be stated as follows: 

Lemma 2 Consider the term obtained by evaluating first the residue at the pole l\ due to 
in the h-integration and then the residue at the pole l^ due to Pij (with 4 inserted) in the 
h- integration. It differs from the one obtained by first calculating the residue due to Pij and then 
to Pi.i (with q inserted) by a sign only. 

To prove it, we write the inverse propagators as Pi : i — on + (3nh +(3i2h and note that the locations 
of poles li , \ 4 and l^ are obtained by solving a linear system 

(£)-(;)♦(£ £)(£W!)- 

Hence the two orders of integration amount to the two ways of solving such a system: first solving 
the first line and inserting it into the second or vice versa. The solutions inserted into the remaining 

(i) (i) (i) (?) 

P q are therefore the same: l\ = l\ and l\ = l\ ■ This, however, specifies only proportionality 
in our second lemma. The proportionality-constant —1 can readily be obtained by writing down 
the residual term after evaluating the (1, i) — > (2, j)-order: 



l ' 3 (/3il/?j2 — Pl.q\,W ,U) ' ' ' ' 

which is antisymmetric in i and j. 

In the planar case, Lemma 2 can always be applied twice to restrict the number of individual 
terms by a factor 4. For the topology HI we thus end up with 9 terms. 

One can further see that the coefficients of the middle variables in Pi are generated by subse- 
quent multiplication of terms linear in the edge variables (the /3y in (Q) are linear in edge variables 
according to (0)) and can therefore be factorized when inserted in the Heaviside functions. There- 
fore the domains of integration in the edge variables are always bounded by a number of linear 
hypersurfaces. 




Figure 2: Constraints after the h- and the ^-integrations 

E.g. when we perform the ^-integration for a scattering-amplitude of type I I I and apply lemma 
1 on all the resulting residues, we obtain the infinite domains sketched on the left of figure ||. (The 
signs denote the weight of the integrand there and stems from the numerator in (||)). After the 
^-integration we obtain the domain in the middle and when multiplied, only a finite triangle in 
the fc3-/3-plane remains. 

Performing the l\- and the fci -integrations in the same way we encounter poles on the real axis. 
They are, however, purely artificial and a consequence of choosing equal imaginary parts rj in all 
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Figure 3: Constraints after the l\- and the ki -integrations in topology lb). The exact parameters 
of the borders are functions of fc 3 and I3 . 



propagators. Treating the integral as a Cauchy Principal Value integral, the residue theorem can 
still be applied in these cases since there are only odd coefficients in the integrand's Laurent-series 
but the residue contributes only with a weight iri instead of 2iri [|l3|| . 

When this procedure is continued and all the constraints are combined applying both lemmata, 
we get the 6 domains of integration in the fco-/o-plane shown in figure || for a scattering-amplitude 
III. (There are only 6 domains remaining instead of 9 due to an incompatibility of constraints 
after the h- and the fc2-integrations in 3 of them.) The exact parameters of the borders of these 
domains depend on k^ and I3. We emphasize again, that the resulting domains of integration in 
the edge variables are always finite and bounded by linear functions — even in the case of crossed 
topologies — thus making them accessible to an unambiguous reparametrization and subsequent 
next integration. In the case of crossed topologies, lemma 2 can be applicable more than twice 
because of the presence of two mixed propagators. The numbers of remaining terms turn out to 
be 12 in I X and 5 in I x . 

3 Towards a 3-dimensional Integral 

The procedure outlined above produces a 4-fold integration of rational function over a finite volume 
in the edge variables: 

r k 3 } r l ° 0) r k ° 0) p(l kn u 

dh dkj di dk 2, 1 , ( 8 ) 

The appearance of a numerator containing integration variables is due to the solutions of P = 
with respect to middle variables inserted into the remaining propagators. After partial-fractioning 
this integrand can always be transformed into one, where Q is no more than quadratic in the 
edge variables. The integration-domains in the fco-Zo-plane which depend on k^ and I3 can be 
mapped into some fixed domain independent of ks and ^3 using a suitable linear transformation. 
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The rational function can then be integrated once more, resulting in other rational functions, 
logarithms and arcustangens where care has to be taken of the small imaginary part and the 
position of branch-cuts. The coefficients of the result will generally involve square roots of quartic 
functions in the remaining three variables. This three-dimensional integration should be accessible 
to numerical evaluation using vegas |l5| , |l6| , or similar routines. 

4 An Example 

We will now test our method on a simple example: a decay-amplitude of type /V where a heavy 
scalar particle decays into a lighter and two massless ones as shown in figure |j. The choice of 
momentum-flow indicated there results in the inverse propagators 

Pis =l 2 -m\ + irj, P l<2 = (l+p 3 ) 2 -mj + ir], 

Pk,3 = (k — pi+ P2) 2 - ml + ir/, P k ,i = (k - pi) 2 -rn\ + ir/, 

Pi+kfi = (l + kf -m\ + ir]. (9) 




Figure 4: Choice of momentum-flow for a decay-process of type III b) 



The most convenient Lorentz frame for this case turns out not to be the rest-frame of the 
decaying particle but the one in which the 3-momenta of the light-like particles are antiparallel 
and aligned to a specific coordinate-axis, say x. This frame can always be reached by a Lorentz 
transformation as long as the two light-like particles do not have equal momenta. 

When integrating out the middle variables with these conditions and applying the rules found 
in section |^ one finds, however, incompatible constraints in the h- and the fci -integrations. To 
get out of the dilemma, the light-cone condition for the two massless external particles must be 
temporarily relaxed]]. A possible choice of momenta is 

Pi = {q V >Px,Py,Q), P2 = (<?7>&»°>°)> PS = (?7> _ fc>°)0)) 

Pa =Pi -P2 -Ps> = (iv -2q 1 ,p x ,p v ,Q), (10) 

where g 7 = lim e ^ + o q x + e is tacitly assumed. 

Now the integrations in the middle variables can be performed resulting in the domain of 
figure § for the fco- and the ^-integrations in addition to the one already known from figure || for 
the k%- and ^-integrations. 

If we use a linear transformation in the ko-lo-p\&ne in order to map this domain into a unit- 
square fco = (0...1), Iq = (0...1), another one to map the triangle in the fc3-/3-plane into a 
unit-triangle ^3 = (0 ... 1), l 3 = (0 . . . 1 — fc 3 ), the limit e — > can safely be performed and we 
obtain the representation 

V( Pl ,P2,p 3 ) = - (2m) 4 [ dk 3 [ dl 3 [ d~k f dlo —~ — ~ l —- (11) 

Jo Jo Jo Jo [pi fco + OoJ to + \a\ fco + ao) + tr) 

1 This difficulty is a manifestation of a degeneracy pointed out in jl|| which prohibits the application of our 
method to two- and three-point functions. 
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where the coefficients are given by 

a = 16 (k 3 l a m\j(l 3 + k 3 - 1) - k 3 ml - l 3 ml 
~ h h (p 2 x +p 2 y -ql + 2 q y (q v - p x )) ) , 

ai = 16 1 3 (2 q 1 k 3 (q n -p x )+ ml -ml), (12) 
b Q = 16 k 3 (2 q 1 1 3 (2 q 1 - q n - p x ) + m\ - ml) , 
h = -64k 3 l 3 q 2 . 

In this case, all terms quadratic in or ko in the denominator have canceled, allowing for an 
integration in both of these variables. Splitting the integral into Principal Value integral and 
^-function separates real- and imaginary part. We obtain a numerically stable two-fold represen- 
tation with dilogarithms in the real part and logarithms in the imaginary part which displays the 
correct behaviour even at thresholds. 

By letting q 7 — > we restrict our kinematical regime even further and obtain a two-point 
sunset-topology with two propagators squared. Numerical stability, however, breaks down as 
we approach this limit. A look at the arguments of the dilogarithms in our two-dimensional 
representation expanded in small g 7 reveals why: 



Urn Re (V( Pl ,p 2 ,p 3 )) = Re (P(p?)) = -(2™) 4 f dk 3 [ dl 3 
" Jo Jo 

Re (Li 2 (l + f 2i qty - Li a (l + r 23 q 2 ,) _ Li 2 (l + f 14 + Li 2 (l + f 13 <%)) . (13) 
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The coefficients in this formula are 

4 k 3 l 3 p\ 



(m| — rn^)[m\ — m 2 ,) 

(k 3 mf + l 3 mfj (1 — k — I) + ml l 3 k 3 
+ 4 = = . (14) 

(m§ ~ ml) (ml - m%)(k 3 + l 3 - 1) 

with pi — q 2 , — pi — Py as the only combination of exterior momenta (c.f. eq. entering P, as 

it should be. 

Numerical stability can be restored in ( fl3"l) if the expansion of the dilogarithms at their branch- 
point as generalized Taylor series 

Re (Li 2 (l + e)) = y + e(l - In |e|) + 0(e 2 ) 
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and the relation f 24 — ^23 — fu + fi3 = are used. 

We compare this with numerical results obtained via another method [Q and find agree- 
ment below and above threshold (figure ||). (For definiteness, all masses have been chosen equal, 
requiring yet another expansion due to the denominators in (14).) 




j I , I , I , I , I , I , I , I , I , I , I , I , L 



1 2 3 4 5 6 7 8 9 10 11 12 13 

p 2 /m 2 



Figure 6: Comparison of our method with a calculation of the scalar sunset-graph with two squared 
propagators of mass m. (Diamonds represent numerical results from our calculation.) 



5 Conclusion 

The method for calculating scalar 2-loop box-functions with arbitrary internal masses proposed 
in Jl[] turns out to deliver a moderate number of 4-dimensional integrals, which can always be 
reduced further to 3-dimensional representations — in some cases even 2-dimensional ones. In 
sample-cases we have been able to produce reasonable numerical results in arbitrary kinematical 
regimes below and above threshold. In the limit of limiting kinematical points, numerical stability 
is lost but can be restored by expanding the representation around that point. We hope to obtain 
similar results for all the 5 genuine 2-loop box-functions and incorporate them into xloops |ft0|| . 
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